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In this paper we study Laurents paper. 

Laurent shows that it is possible to construct Phase modulation by the superposition of amplitude modulated pulses. 
The new idea : 

We can separate the method of the superposition from the shape of the amplitude modulated pulses. We can then 
use different pulses of the same duration and /or fewer of them to construct modulation schemes that have desirable 
properties in the spectral and approximate to constant envelope schemes. 

In particular we can use pulses with narrower bandwidth (but of the same support) to generate a modulation scheme that 
has a narrower bandwidth with a corresponding increase in the variation of the envelope of the signal 
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In this paper we effectly study Laurents paper. 

We develop the modulating function that is used in GSM first to use as an example to demonstrate Laurent's idea. 



3 

T : = 



812500 
BT : = 0 . 3 



VLog[2] 



n[t_] := 



2 7T BT 



f 2 7T a T 



Ideally we would and did define the effect of the convolution of h[t] by the formula below. However, this version of 
Mathematica gives an error. 

T 

kfii [t_] := Release [Module [{r>, j~ T h I fc ~ r ? ^ 

We get round this by numerically fitting an interpolation function. 

kfii [PulseWidth_] [t_] := Module [{xl ; x2, x3}, xl = 



(N[#l, 40]&) [Table[{tl, h[tl " r ? dr } # r fcl 

x2 = Interpolation [xl] ; x2 [t] ] 
Plot[h[tT], {t, -5, 5}, PlotRange -> All] 



PulseWidtn PulseWidth T 



20 



}]] 



- Graphics - 

In this section we define (f) L(t) the phase modulation function with a pulse width of length L at time t 

Modulationlndex := — 
2 

We use the definitions in the paper 

« : = N [Modulationlndex tt] 

C : = Cos [£] ; 
S := Sin[ff] ; 
J := //Cliop); 
M := 2 L ' 1 ; 

Clear [PhaseAngle] 

PhaseAngle [L_] [t_] / ; t <; 0 : = 0 

PhaseAngle [L_] [t_] /;t*X«T:=ff 
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PhaseAngle[L_] [t_] : = 

PhaseAngle[L] [t_] = Module [{xl, x2, x3, x4, x5, x6}, xl=h £il [3LT] f tl - — 1; 

2 J 

x2 = Table[{t2, Bj Evaluate [xl] dltl} , {t2 # 0, LT, -^"}] ; Interpolation [x2] [t J ] 

RuleDelayed: : rhs : Pattern t_ appears on the right-hand side of 

rule 0 L _,t_ » [<t>L.z_ = Module [{xl, x2, x3, x4, x5, x6}, xl = h fil [3 LT] [tl - ; 
x2 = Table [ «1» ] ; Interpolation [x2 ] [t] ] J . 

When running the notebook we fix the value of L here and the rest of the graphs and functions use this value 
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L := 8; 

PhaseFunction[t_] = 0 L ,t/ 

phasepoints = Table [{t, PhasePunction[t T] } // N, {t, 0, L , 1/40}]; 
Lis t Plot [phasepoints, Plot Joined -> True] 

NIntegrate: :ncvb : NIntegrate failed to converge to prescribed accuracy after 7 
recursive bisections in tl near tl = -3 . 68196 x 10~ 7 . 

NIntegrate: :ncvb : NIntegrate failed to converge to prescribed accuracy after 7 
recursive bisections in ,tl near tl = -7 .4914 x 10" 7 . 

NIntegrate: :ncvb : NIntegrate failed to converge to prescribed accuracy after 7 
recursive bisections in tl near tl = - 5 . 43232 x 10" 7 . 

General: : stop : 

Further output of NIntegrate: :ncvb will' be suppressed during this calculation. 

NIntegrate: :slwcon : Numerical integration converging too slowly; suspect singularity, value 
of the integration is 0, oscillatory integrand, or insufficient WorkingPrecision. 
If your integrand is oscillatory try using the option Method->Oscillatory in NIntegrate. 

NIntegrate: : slwcon : Numerical integration converging too slowly; suspect singularity, value 
of the integration is 0, oscillatory integrand, or insufficient WorkingPrecision 
If your integrand is oscillatory try using the option Method->Oscillatory in NIntegrate. 

NIntegrate: : slwcon : Numerical integration converging too slowly; suspect singularity value 
of the integration is 0, oscillatory integrand, or insufficient WorkingPrecision 
■If your integrand is oscillatory try using the option Method- >Oscillatory in NIntegrate. 

General :: stop : 

Further output of NIntegrate: : slwcon will be suppressed during this calculation. 
1.5; 

1.25 • 
1 ■ 

0.75- 

0.5 
0.25 



Graphics 



We can define an phase modulated signal to be 



As can be seen from the graph 0 t is 0 for most of the time and is nearly 4> after LT . We denote by <£ L , t the pulse that is 
set to these values 
i.e 0 L tt =0 t<0; 

* L ,t =«>t >LT; 

0L,t = 0t Otherwise 



S t = E 1 ^— *" n ^-nT 

We now consider t = N T + r where 0 < r < T 

S t = E I ^"=-~ ffn ^ NT+T -° T 
S t = E 1 ff " £l.t + (N-o>T 



Making use of the fact that <f> hi t = 0 t < 0 ; 
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S t = E 1 ^n=-« *U r + (N-n) T SjL N-L + 1 a » *Ur+ (N-n)T 

Making use of 0 L/ 1 = 4> t > L T we get 

S t = E 11 ---^ * E I2: n=N-L+l ff n *L,r + (N-n>T 

And using the notation for J : = * v ' 
Making the substitution n = N - i 
Express this as 

S t = J ^ 1 (El«N-l*L.r + JT) 

i = 0 

Since we have 

( Sin[A - B] Sin[B] ( - ) > 

1 Sin[A] + * il^f Sia[X " ] II J " Si "*> lif * 

E -I B 

And 
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(Sin[A-B] Sin[B] <e* x -e-J*)\ 

1 Sin [A] * * "iZnUT SintX " 3 Jj J " Sia *> lif * 

We can make the following substitutions S := Sin[4>] 

v s s J 

We can make the substitution 

t_] /; 0< t < LT := 0 L , t 

t_] /;2LT>tiLT := ff — t - l t 
We need to put this to avoid the function being extrapolated 

*[ L - # t_] /; ! (0<t<LT) &&! (2 LT > t 2s LT) : = 0 
Clear [torn, xl] 

torn = (xl[t_] = *[L, tT]; Table [{t, xl[t ]} / / N, {t, 0, 2L, 1/40}]); 
Lis tPlot [torn. Plot Joined -> True] 

Now define 

Clear [Laurent S] 

Laurent S [L_] [n_] [t_] := Sin[ #[L, t + nT]] / S 
Plot [LaurentS[L] [2] [t T] , {t, -2, 10}] 

- Graphics - 

* = *> I S S J 



*=o I S S J 



s t = L n 1 f Sin[ ^ t - NT+iT ^ LTl . T a N ., Sin ^t-NT. iT ] ^ 

«=o I S S J 

S t = J L n \ LaurentSp + L - N][t] + J ff "-' LaurentS[i - N][t] ) 

We try to get an insight into n ( LaurentSfi + L - N][t] + J*"-' LaurentSfi - N][t] ) 
To this end we try defining 

Temp[L_, N_] [t_] : = J~J (5 [i + L - N] [t ] + J*"-* 5[i-N] [t] ) 



v 
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(Temp [4, 0] [t] // Expand // Chop) /. Complex[0, 1.] -> Complex [0 , 1] 

ia T a :r a + V a ° Sl ° ] [t] Sll] [t] Sl2] [t] ^ [3] [t] ^I 0 - 3 ^ 2 ^ 1 S[l] [t] S[2] [t] S[3] [t] 5[4] [t] + 

C 3 ;r 5[o] [tj 5 ^ sis] [t] + i— > [t] s[3] t ^ 4 1 5 6 j + 

5[0] [t] s[2] [t] 5 [5] [t] 5[6] [t] + l°-» 5[3] [t] S[4] [t] 5[5] t] 5 6] t] V 
• ^^.co 5[0] [t] 3[1] [t] 5[2] [t] [t] + [t] s[2] s L {. ' 

X a- 2 .a 0 5[Q] [fc] 5[2] [fc j ^ [5] [t] s[1] [t] + ia . 2 5[2] [fc 5[ l + J 

ja-i^ao 5 [0] [t] S[l] [t] S[6] [t] 5[7] [t] + [t] 5[4] [t] 5[6] [t] 5[7] [t] + 

I a ° 5[0] [t] S[5] [t] S[SH-t] S[l] [t] + 5[4] [t] S[5) [t] S[S) [t] S[7) [t] 

In the next section we define the 

Notation [a^^i^ » AlphaKI [LL_] [K_ # ii_] ] 

Clear [AlphaKI] 

AlphaKI [LL_] [K_, i_] /; 0<i<LL&&0^K< 2"" 1 : = 
Module [{xl, x2, x3, KNum} , xl : = ^x2 = Mod [KNum, 2]; KNum = * Num " x2 ; x2 j ; 
KNum = K; x3 = Table [xl, 0, LL - 1}] ; x3 ] 

Check the function. 

Sum[2*(ii - 1) AlphaKI [6] [13, ii], {ii, 1, 6 -1}] 
13 

We can now define the LaurentC functions. 
Clear [LaurentC] 

General: rspelll : Possible spelling error: new symbol name "LaurentC- is similar 
to existing symbol "LaurentS". 

LaurentC [L_] [K_] [t„] / ; 0 i K < 2 L : = 

L-l 

LaurentS [L] [0] [t] [~] LaurentS [L] [ii + L AlphaKI [L] [K, ii] ] [t] 

ii=l 

LaurentLK gives the support of the C functions 

LaurentLK [L_] [K_] ' = 
Module [{xl}, xl = Table [L (2 - AlphaKI [L] [K, ii] ) - ii, {ii, l, L-l}]; Min[xl]] 

Table [LaurentLK [L] [ii], {ii, 0, 2 L - X -1}] 

{9, 7, 6, 6, 5, 5, 5, 5, 4 # 4, 4, 4, 4, 4, 4, 4, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 

3, 3, 3, 3, 3, -2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2 

2 2 2 2 2 2 2, 2, 2, 2, 1, 1, l, 1, 1, 1, 1, i, 1, lf i, i 2 i x x {' 

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, i,.i, i, i, i/i, 1§ i i ! i x 
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, X/ l, i, 1, 1} A ' ±9 ±f 

CFunctionPlot [K_] := Module [ {xl, x2 }, xl [t_] = LaurentC [L] [K] [t T] ; 

x2 = Table [{t, xl[t]}, {t, 0, LaurentLK [L] [K] , 1/40}]; ListPlot[x2, 
PlotJoined -> True, PlotLabel -> StringJoin["\n \n C = ", ToString[ii] ] 
FormatType -> StandardFora, PlotRange ->All]] 

Plot at most the first 7 C functions 



Table [CFunctionPlot [ii] , {ii, 0, Min[6, 2 L " 1 -1]}] 



0.8 
0.6 
0.4 
0.2 



0. 07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 



0. 0014 
0.0012 
0.001 
0.0008 
0.0006 
0.0004 
0.0002 



0.0014 
0.0012 
0.001 
0.0008 
0.0006 
0.0004 
0.0002 
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2 .5xl0" 6 
2xl0~ 6 

l.SxlO' 6 
lxlCf 6 
5xl0" 7 



7xl0" 6 
6xl0" 6 
5x10" 6 
4xl0" 6 
3xl0" 6 
2xl0" 6 
lxlO" 6 



5xl0" 7 
4xl0" 7 
3xl(T 7 
2xl0" 7 
lxlO" 7 
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{- Graphics - , - Graphics - , - Graphics - , - Graphics - , - Graphics - , - Graph 
- Graphics - } 

In this section we define the modulating function. We intend to build the modulator from first principle: 

y 1 ^ 0 **' 0 c K [t-nT] 
~T — 'k=o 

n--oo 

Change the order of summation 

M -1 9 

Given that we are interested at t = N T + t 
m-1 a 

V J C K [r + (N - n) T] 

Change variable n n' + N 

M -1 m 

Notice that the support of C K is positive only. 

M -1 0 

Using LaurentLK as 

^^f 1 n' = -Laureate [L] [K] + X 

Changing the sign of n 

M _1 [ LaurpptUt [L»] [K] - 1 



Let the value of A 0 , o be specified. 

Then we can calculate the value of A 0 , _i ,A 0 , _ 2 ,A 0 , -3 using the following algorithm 

^n-l = Ao, n ~ On 
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BitSeq = Table [1, {i, 1, 20} ] 

{1, 1, 1, 1, 1, 1, 1, 1, 1, 1, l, i, i, i, i, i, x, lf lf 1} . 

ANKInitialStateSetUp[L_] [K_] [InitBitSeq_, AccumulatedPhase^] : = 
Module [{xl, x2, x3, x4, x5, aeuphase, initbitseq}, 
initbitseq = InitBitSeq; 
aeuphase = AccumulatedPhase; 
UpdateSeq : = 

Module[{}, xl = aeuphase - Sum [initbitseq [ [i ] ] AlphaKI [L] [K, i] , {i, 1, l - m - 
aeuphase = aeuphase - First [initbitseq] ; initbitseq = Rest [ initbitseq]; xl]; 
Table[ UpdateSeq, {i, 1, LaurentLK[L] [K]}]] 

With the BitSeq we have we expect a decreasing sequence starting at 0. 
ANKInitialStateSetUp[L] [0] [BitSeq, 0] 

{0, -l,.-2, -3, -4, -5, -6, -7, -8} 

Testlnit = Table [1, {i, 1, l}] 

{1, 1, 1, 1, 1, 1, 1, 1} 

AKN[L_i [K_] [{State., AccuarulatedPhaseJ] : = 
AccumulatedPhase - Sum [State [ [i + 1] ] AlphaKI [L] [K, i] , {i, 1, L -1}] 

Options [Modulator] := {StartingQuadrant 0, InitialState Table [1, {i 1 20U 
Samplinglnterval - T/32, NumberOf Curves 4, ModulatingPulse Laure'ntC} 

BitSeq = Table [1, {i, 1, 20}] 

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, i, l, i, i, i, i, i, 1} 

Modulator [L__] [BitSeq..., Opts ] : = 

Module [ 

{xl, x2, x3, x4, x5, x6, state, AccumalatedPhase, seq, AKNState, Curves, Pulse}, 
xl = Samplmglnterval/. {Opts} /. Options [Modulator] ; 
state = InitialState /. {Opts} /. Options [Modulator] ; 
x3 = StartingQuadrant/. {Opts} /. Options [Modulator] ; 
x4 = Samplinglnterval/. {Opts} /. Options [Modulator] ; 
Pulse = ModulatingPulse /. {Opts} /. Options [Modulator] ; 
Curves = (NumberOf Curves / . {Opts} /. Options [Modulator] ) -1; 
AccumulatedPhase = x3 - BitSeq[[l]]; 
seq = BitSeq; 
Table [AKNState [K] = 

ANKInitialStateSetUp[I.] [K] [state, AccumulatedPhase], {K, 0, Curves}! • 
x5 := Module [{}, state = Join [ {First [seq] } , Drop [state, -1]]; 
AccumulatedPhase = AccumulatedPhase + First [seq]; 
seq = Rest [seq]; 

Table [AKNState [K] = Join [ {AKN [I,] [K] [ {state, AccumulatedPhase} ]} , 

Drop [AKNState [K] , -1]], {K, 0, Curves}]; 
X6[r_] = Sum[Sum[ (J) estate [KJ [t i + m pulse [L] [K] [r + i T] , 

{i, 0, LaurentLK[L] [K] -1}], {K, 0, Curves}]; 
Table[x6[r], {r, 0, T - x4, x4}]]; 
Table [x5, {kk , 1, Length [BitSeq] } ] //Flatten] 

General: : spelll : 

Possible spelling error: new symbol name "state" is similar to existing symbol, - State" . 

General: : spelll : Possible spelling error: new symbol name "AccumalatedPhase" is 
similar to existing symbol "AccumulatedPhase". 
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RandomBitSeq = Table [Random [Integer, {0, 1} ] # {i, 1, 40)] //Map[# (-2) + 1&, #]& 
Modulator [I*] [BitSeq, Samplinglnterval T/10]; 

Modulator [L] [BitSeq, Samplinglnterval T/10, NumberOf Curves -> 8]; 
Modulator [I*] [BitSeq, Samplinglnterval -> T/10, NumberOf Curves -> 1]; 
torn = Modulator [L] [RandomBitSeq, Samplinglnterval T/10, NumberOf Curves -> 2] ; 
Lis tPlot [Im[ torn] , Plot Joined -> True] 

1 



-0.5 




- Graphics - 



•c 



LaurentPaperPatents. nb 



12 



I*i st Plot [Abs [torn] , Plot Joined -> True] 



1.002 
1.001 



0.999 
0.998 



m 



23 



M 



M 



DC 



w 



400 



- Graphics - 

{Re [torn], Im[tom]> // Transpose // ListPlot [#, Plot Joined -> True, AspectRatio -> 1]& 




- Graphics - 

torn = Modulator [L] [RandomBitSeq, Samplinglnterval -+ T/10, NumberOf Curves -> 1] 
Li s t Plot [Im[ torn] , Plot Joined -> True] 




400 



Graphics - 
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List Plot [Re [torn] , Plot Joined -> True] 




400 



- Graphics - 

ListPlot [Abs [torn] , Plot Joined -> True] 



1.06 
1.04 
1.02 



0.98 
0.96 
0.94 



100 



200 



M 



00 



- Graphics - 
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{Re[t©m], Im[tom3> // Transpose // ListPlot[#, PlotJoined -> True, AspectRatio -> 1]& 









fz 

( 0 • 2 5 




> 


-0.75-0.5 -0.25 
i -0.25 


0.25 0.5 0.75 
^^^^ 







- Graphics - 
Calculating the bandwidth of the signal 

xll = Table [LaurentC [L] [0] [t TJ , {t, 0, L + 1 - 1/200, 1/200}]; 
((Drop [xll, -1] - Drop [xll, 1]) 200 )*2 / 200 // Apply [Plus, #]& 

1.29684 

Here the pulse 6T is just stretched to 8T to see the performance 



LaurentPaperPatents.nb 



PhaseFunction[t_] = # 6 , t ; 

phasepoints = Table [{t, PhaseFunction[t T] } // N, {t, 0, 6, 1/40}]; 
ListPlot [phasepoints. Plot Joined -> True] 

NIntegrate: : ncvb : NIntegrate failed to converge to prescribed accuracy after 7 
recursive bisections in tl near tl = - 1 . 91796 x 10~ 7 . 

NIntegrate: : ncvb : NIntegrate failed to converge to prescribed accuracy after 7 
recursive bisections in tl near tl = - 5 . 52869 x 10" 7 . 

NIntegrate: islwcon : Numerical integration converging too slowly; suspect singularity value 
of the integration is 0, oscillatory integrand, or insufficient WorkingPrecision 
If your integrand is oscillatory try using the option Method- >Oscillatory in NIntegrate. 

NIntegrate: : ncvb : NIntegrate failed to converge to prescribed accuracy after 7 
recursive bisections in tl near tl = - 1 . 9669 6 x 10" 7 . 

General : : stop : 

Further output of NIntegrate: : ncvb will be suppressed during this calculation. 

NIntegrate: islwcon: Numerical integration converging too slowly; suspect singularity value 
of the integration is 0, oscillatory integrand, or insufficient WorkingPrecision 
If your integrand is oscillatory try using the option Method->Oscillatory in NIntegrate 



Llatory in NIntegrate. 
General: : stop : 

Further output of NIntegrate: : slwcon will be suppressed during this calculation. 

- Graphics - 

TestPulse[8] [0] [t_] = LaurentC [ 6 ] [0] [t 7 / 9 ] 

^Aco^-Y-]]^Ms. i^JlsM^, 812T00" + ^T~]] 
sin W 6 ' 20^25-* ^]] sin W 6 - leAoo ♦■¥-]]*!»[*[«. 2*.]] 

Plot[TestPulse[8] [0] [t] , {t, 0, 9T}, PlotRange -> All] 

- Graphics - 

RandoxoBitSeq = Table [Random [Integer, {0, 1}], {i, 1, 100}] //Map[# (-2) + 1&, #]& 
{1, 1, -1, -1, -1, 1, 1, -1, 1, -1, 1, -l, 1, -l, -i, -i, ! ! 1# * n .iiii 

i, -l, i, -l, i, -l, i, i, -i, -i, -i, i, -i, i # i, lf lt j lt j lt '_i i ti t ±9 
1, i, i, -l, i # -l, -i, i, i, i # -i, -i, i, i, i, i, i, i, 1# x§ lf _ lt _ lr lf : Xt lt _ 1} 

torn = Modulator [L] [RandomBitSeq, Samplinglaterval -> T/10, NumberOf Curves -> 1 
ModulatingPulse -+ TestPulse] ; ' 

ListPlot [Im[ torn] , Plot Joined -> True] 

- Graphics - 

ListPlot [Re [torn] , Plot Joined -> True] 



- Graphics - 



LaurentPaperPatents.nb 

ListPlot [Abs [torn] , Plot Joined -> True] 

- Graphics - 

{Re [torn], Im[tom]} // Transpose // ListPlot [#, PlotJoined -> True, AspectRatio -> 1]& 

- Graphics - 

In this section we consider filtering the pulse 
Needs ["bessel*"] 

fil = BesselFilter[8] [LowPass3dB [1 / T 2 Pi 0 . 3] ] [Frequency Response] [s] ; 

Plot [20 Log[10, Evaluate [ (fil /. s 2 Pi I f / T) //Abs]], {f, 0, 2}, 
GridLines -> Automatic, Axes Label -> {"1/T", "dB"}] 

- Graphics - 

Res[t_] = InverseLaplaceTransform[fil, s, t] 

2027025 ( (9.36812 + 0. 1) e" 897172 ' * Cos [139301. t] - 

(13 .7944 + 0. I) e~ 835672 • t C os[420044. t] + (4.73111 + 0. I) E" 701358 * c Cos[708768. t] - 
(0.3 04883 + 0.1) E" 45581S - 6 Cos [ 1 . 02 016 x 10 6 t] + 
(29.1264 + 0. I) E" 897172 * z Sin[139301. t] - 

(9.89398 + 0\ I) E" 835672 - * Sin[420044. t] - (0.321205 + 0. I) E' 701358 - * Sin[708768. t] + 
(0.375154 + 0. I) E" 455818 * t Sin [ 1 . 02 016 x 10 6 t] ) 

Plot [Res [t T] , {t, 0, €}, PlotRange -> All] 

- Graphics - 

filTaps = Table [Res [t T] , {t, 0, 6, 1/20}] // Chop; 

We estimate the group delay to be 1 .8T 
Length. [ filTaps ] 
121 

pulse = Table [LaurentC [8] [0 ] [tT], {t, 0, 9, 1/20}]; 

General: :spelll : 

Possible spelling error: new symbol name "pulse" is similar to existing symbol "Pulse". 
pulse2 = Table [LaurentC [8] [1] [t T] , {t, 0, 7, 1/20}]; 
Needs [ " Convolve * " ] 
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FilteredPulse = T/20 MapConvolve [pulse, filTaps] //Take[#, {34, Length[#] }]& 

. {1.73731X 10" 15 , 5.35887xl0" 15 , 1 . 60413 x 10" 14 , 4 . 6627 8 x lCT 14 , 1 . 31688 x 10" 13 , 
3.61578 xl0" 13 , 9.65729xl0" 13 , 2 . 51044 x 1CT 12 , 6 . 3 55 07 x 10" 12 , 1 . 56747 x 10" 11 , 
3.76889 X 10-- 11 , 8 . 83867 x 10" 11 , 2 . 02275 x 10~ 10 , 4 . 51958 x 10' 10 , 9 . 86455 x 10" 10 , 
2 .10426 x 10" 9 , 4.38922 xl0* 9 , 8 . 95697 x 10 -9 , 1 . 78914 x 10~ 8 , 3.4999xl0 -8 , 6 . 70842 x 10 -8 , 
1.26054xl0" 7 , 2.32319xl0' 7 , 4 . 20166 x 10~ 7 , 7 . 46062 x 10~ 7 , 1 . 30124 x 10" 6 , 
2.23035xl0~ 6 , 3.75853 xl0" 6 , 6.23xl0 -6 , 0.0000101618, 0 . 0000163174, 0 . 0000258049 
0.0000402068, 0.000061746, 0.000093496, 0.00013964, 0.000205786, 0 000299337 
0.000429924, 0.000609895, 0.000854853, 0.00118424, 0.00162196, 0 00219697 
0.0029439, 0.00390362, 0.00512368, 0.00665873, 0.00857068, 0.0109288, 0 0138093 
0.0172954, 0.021476, 0.0264451, 0.0323001, 0.0391407, 0.0470666, 0.0561754 
0.0665603, 0.0783079, 0.0914948, 0.106186, 0.12243, 0.140261, 0.159691, 0.180711 
0.20329, 0.22737, 0.25287, 0.279682, 0.307673, 0.336686, 0.366539, 0 397029 

. 0.427934, 0.459013, 0.490013, 0.520668, 0.550708, 0.579856, 0.607837, 0.634379 
0.65922, 0.682107, 0.702805, 0.721098, 0.73679, 0.749713, 0.759727, 0.76672 
0.770616, 0.771369, 0.768969, 0.763442, 0.754846, 0.743274, 0.728852, 0 711736 
0.692108, 0.670178, 0.646176, 0.62035, 0.592964, 0.56429, 0.534607, 0.504193, ' 
0.473326, 0.442277, 0.411305, 0.380656, 0.350558, 0.321221, 0.29283, 0.26555 
0.239519, 0.214852, 0.191636, 0.169936, 0.149791, 0.13122, 0.114217, 0.0987613, 
0.0848118, 0.0723143, 0.0612016, 0.0513967, 0.0428146, 0.035365, 0.0289542 
0.023487, 0.0188687, 0.0150065, 0.0118109, 0.00919671, 0.00708412, 0.0053992 
0.00407439, 0.00304877, 0.00226811, 0.0016848, 0.00125761, 0.000951417, 0 000736739 
0.000589309, 0.000489543, 0.000422011, 0.000374899, 0.000339478, 0.0003096 
0.000281223, 0.000251974, 0.000220764, 0.000187441, 0.000152497, 0.000116824 
0.0000815146, 0.0000477044, 0.0000164566, -0.000011319 8, -0.0000349206 
-0.0000538729, -0.00006793 91, -0.0000771062, -0.0000815622, -0.0000816651 
-0.0000779073, -0.000070877, -0.000061222, -0.0000496142, -0.00003 67186, 
-0.0000231668, - 9 . 53584 x 10~ 6 , 3 . 66817 x 10" € , 0.0000160205, 0.000027183 8, 
0.0000369093, 0.0000450345, 0.0000514769, 0.0000562268, 0.00005933 63, 0.000060909 
0.0000610882, 0.0000600456, 0.0000579702, 0.0000550589, 0.0000515076, 0 0000475036 
0.0000432203, 0.000038813, 0.0000344153, 0.0000301387, 0.0000260716, 
0.0000222801, 0.0000188096, 0.000015687, 0.0000129228, 0.000010514, 
8.4469 xl0" 6 , 6.69967 xl0" 6 , 5 . 24481 x 10' 6 , 4 . 05145 x 10 -6 , 3 . 08723 x 10' 6 , 
2.3199xl0' 6 , 1.71859 xlO" 6 , 1 . 25466 x 10" 6 , 9.0236xl0" 7 , 6 . 39109 x 10" 7 , 
4.45602 xl0* 7 , 3.05722 xl0" 7 , 2 . 063 18 x 10 -7 , 1 . 3 6898 x 10" 7 , 8 . 92721 x 10~ 8 , 
5.71866xl0 -8 , 3.59686xl0" 8 , 2 . 22017 x 10" 8 , 1 . 34418 x 10" 8 , 7 . 97 812 x 10" 9 , 
4.63953 xl0" 9 , 2 . 64191 x 10" 9 , 1.4722 xl0" 9 , 8 . 02333 x 10" 10 , 4 . 273 88 x 10" 10 , 
2.2234X10" 10 , 1.12886 xl0- 10 , 5 . 5927 6 x 10' 11 , 2 . 70717 x 10" 11 , 1 . 28211 x 10' 11 , 
5.92683 xlO" 12 , 2 . 65259 x 10" 12 , 1 . 14156 x 10" 12 , 4 . 73562 x 10" 13 , 1 . 94329 x 10" 13 , 
7.95416xl0" 14 , 3.14517xl0" 14 , 1 . 23089 x 10" 14 , 4 . 73302 x 10" 15 , 1 . 66595 x 10' 15 , 
4.65283xl0' 16 , 1 . 93314 x 10" 16 , 9 . 8773 8 x 10" 17 , 2 . 3 1472 x 10' 18 , - 1 . 50823 x 10" 17 , 
-1-01228x10"", -6.49877xl0~ 18 , - 2 . 08455 x 10" 18 , 1 . 80177 x 10' 19 , 1 . 70062 x 10" 19 , - 
1.9917xl0" 19 , -1. 14008x 10-", - 1 . 66045 x 10" 20 , 3 . 7877 x 10~ 21 , 1 . 47752 x 10" 21 , 
1.94753 xlO" 21 , 3.72298xl0" 22 , 2 . 02535 x 10" 22 , 1 . 36007 x 10" 22 , 6 . 72514 x 10" 24 , 
7.07296xl0" 23 , 6.84862xl0' 23 , 3 . 68916 x 10' 23 , 2 . 6013 x 10' 24 , - 1 . 512 51 x 10" 24 , 
-5.51867 x 10' 25 , 7.52915 xl0" 26 , 5 . 81584 x 10 -25 , 1 . 75371 x 10" 25 , 1 . 22392 x 10' 26 , 0} 
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FilteredPulse2 = T/20 MapConvolve [pulse2, filTaps] //Take[#, {50, Length[#] }]& 

. {1.27972xl0~ 7 , 2.18306 X 10" 7 , 3 . 6 6972 x 1(T 7 , 6.0806 X 1CT 7 , 9 . 93427 x 10" 7 , 1 . 6007 6 x 10" 6 , 
2.54471xl0" 6 , 3.99199xl0" 6 , 6 . 18149 x 10' 6 , 9 . 45059 x 10" 6 , 0.0000142 688, 0.0000212804 
0.0000313562, 0.000045 6572, 0.0000657079, 0.0000934812, 0.000131493 0 000182902 
0.000251616, 0.000342391, 0.000460925, 0.000613932, 0.000809185, 0 00105553 
0.00136282, 0.00174186, 0.00220418/0.00276183, 0.00342705, 0.00421186, 0 00512763 
0.00618453, 0.00739101, 0.00875319, 0.0102744, 0.0119544, 0.0137895, 0 0157714 
0.0178876, 0.0201213, 0.0224509, 0.024851, 0.0272921, 0.0297419, 0.0321652 
0.0345256, 0.0367857, 0.0389085, 0.0408581, 0.0426011, 0.044107, 0 0453496* 
0.0463073, 0. 0469636, *0".""0473081 , 0.0473359, 0.0470483, 0.0464525, 0 0455611 
0.0443919, 0.0429671, 0.0413131, 0.039459, 0.0374364, 0.0352785, 0.0330189,' 
0.0306913, 0.0283283, 0.0259611, 0.0236189, 0.021328, 0.019112, 0.0169912 
0.0149823, 0.0130987, 0.0113504, 0.0097438, 0.00828249, 0.00696697, 0.00579518 
0.00476282, 0.00386365, 0.0030899, 0.00243259, 0.00188191, 0.00142752, 0 00105885 
0.000765388, 0.000536896, 0.000363597, 0.000236336, 0.000146695, 0.0000870728 
0.0000507313, 0.0000318118, 0.0000253258, 0.0000271231, 0.000033 8422, 0.0000428466 
0.0000521516, 0.0000603451, 0.0000665053, 0.0000701193, 0.0000710042, 0 0000692325* 
0.0000650647, 0.0000588885, 0.00005116 65, 0.0000423925, 0.000033 0554, 0.0000236122' 
0.0000144671, 5.95924 xl0' 6 , - 1 . 64542 x 10" 6 , - 8 . 15591 x 10" 6 , -0.0000134543 
-0.0000174907, -0.0000202751, -0.0000218687, -0.0000223734, -0.0000219217 
-0.000020666, -0.0000187698, -0.0000163988, -0.0000137138, -0.0000108643,' 
-7.98457 xl0 -6 , - 5 . 18972 x 10" 6 , - 2 . 57416 x 10" 6 , - 2 . 1065 x 10' 7 , 1 . 84943 x 10" 6 ', 
3.57494xl0" 6 , 4 . 953 19 x 10" 6 , 5 . 98756 x 10" 6 , 6 . 69487 x 10" 6 , 7 . 10243 x 10" 6 , 
7.24522 xlO" 6 , 7 . 16306 x 10' 6 , 6 . 8 9804 x 10" 6 , 6 . 49223 x 10" 6 , 5 . 98586 x 10" 6 , 
5.41573xl0~ 6 , 4.81423 xl0~ 6 , 4 . 20863 x 10" 6 , 3. 6208 xl0' 6 , 3 . 06728 x 10' 6 , 2 . 55958 x 10" 6 , 
2.10471xl0" 6 , 1.70583 X 10' 6 , 1 . 3 62 97 x 10" € , 1 . 07376 x 10" 6 , 8 . 3 4162 x 10" 7 , 6.3908 xl0- ? i 
4.82884xl0' 7 , 3 . 59854 x 10~ 7 , 2 . 64486 x 10" 7 , 1 . 91715 x 10" 7 , 1 . 37044 x 10" 7 , 
9.65984xl0' 8 , 6 . 71324 x 10" 8 , 4 . 59913 x 10" 8 , 3 . 10539 x 10" 8 , 2 . 06612 x 10~ 8 , 1.3542 xl0" 8 , 
8.74155xl0 -9 , 5.55593 xl0" 9 , 3 . 47564 x 10" 9 , 2 . 13 956 x 10" 9 , 1 . 2 9613 x 10' 9 , 
7.73129X10' 10 , 4.54101X10" 10 , 2 . 62031 x 10" 10 , 1 . 478 88 x 10" 10 , 8 . 15113 x 10" 11 , 
4.40884X 10- 11 , 2.37341X10" 11 , 1 . 2 6865 x 10' 11 , 6 . 5895 8 x 10" 12 , 3 . 3 4992 x 10~ 12 ,' 
1.6239 xl0' 12 , 7.01037xl0" 13 , 2 . 33 03 6 x 10" 13 , 7 . 7505 1 x 10" 14 , - 3 . 22249 x 10' 15 ', 
-8.12875 x10-", -7.89371xl0" 14 , -5 . 29088 x 10" 14 , - 3 . 142 97 x 10" 14 , - 8 . 2 1475 x 10' 15 , 
4.12557xl0 -15 , 1.9906 xl0" 15 , 1 . 82 06 8 x 10" 15 , - 4 . 3 6969 x 10" 15 , 4 . 67143 x 10' 16 , 
7.65162 X10" 16 , 5.67707 xlO" 16 , 4 . 933 42 x lO' 16 , 7 . 52562 x 10" 17 , 7 . 756 06 x 10' 17 ,' 
4.50401xl0" 17 , -2.18469xl0" 17 , 1 . 39012 x 10" 16 , 1 . 92 1 64 x 10" 16 , 6 . 31946 x 10" 17 , 
-6.72328xl0 -18 , 1 . 24704 x 10" 17 , 1 . 42 063 x 10' 17 , 7 . 33335 x 10' 18 , 1 . 13717 x 10" 19 \ 
3.90929 xl0 -19 , 1.30285 xl0" 19 , 0} 

FiltPulse[8] [0] = {(Table[tT, {t, 0, 20, 1/20}] // N) // 

Take[#, Length [FilteredPulse] ] &, FilteredPulse} //Transpose// 
Interpolation 

InterpolatingFunction[ { {0, 0.0000492923}}, <>] 

FiltPulse[8] [1] = { (Table [t T, {t, 0, 20, 1/20}] // N) // 

Take[#, Length [FilteredPulse2 ]] &, FilteredPulse 2} //Transpose// 
Interpolation 

InterpolatingFunction[ { {0, 0.0000389538}}, <>] 
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Plot [Filt Pulse [8] [6] [tT], {t, 0, 9}, PlotRange -> All] 




2 4 6 8 

- Graphics - 

Plot[LaurentC[8] [0] [tT], {t, 0, 9}, PlotRange -> All] 




2 4 6 8 

- Graphics - 
Show[%, %%] 




2 4 6 8 

- Graphics - 
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Plot[FiltPulse[8] [1] [t T] # {t, 0, 7},, PlotRange -> All] 



0.04 
0.03 
0.02 
0.01 




1 2 3 4 5 6 7 
- Graphics - 

Plot[LaurentC[8] [1] [t T], {t, 0, 7}, PlotRange -> All] 

0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 




- Graphics - 

Bandwidth [TestPulse_, Supp_] s = 
Module [{xll}, xll = Table [TestPulse [t T] , (t, 0, Supp - 1/200, 1/200}]; 
((Drop[xll, -1] -Drop[xll, 1]) 200 ).* 2 / 200 // Apply [Plus, #]&] 

Bandwidth [FiltPulse [8] [0] , 9] 

0.711659 

BandWidth[I,aurentC[8] [0], 9] 

1.29684 

Bandwidth [FiltPulse [8] [1] , 7] 

0.00338446 

torn = Modulator [L] [Table [1, {i, 1, 100}], Samplinglnterval T/10, 
NumberOfCurves -> 2, ModulatingPulse FiltPulse]; 

torn = Modulator [L] [RandomBitSeq, Samplinglnterval -* T/10, NumberOfCurves 
ModulatingPulse FiltPulse] ; 

ListPlot [Im[tom] , Plottfoined -> True] 



- Graphics - 
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ListPlot [Re [torn] , Plot Joined -> True] 

- Graphics - 

ListPlot [Abs [torn] , PlotJoined -> True] 

- Graphics - 

{Re [torn], Im[tom]} // Transpose // ListPlot [#, PlotJoined -> True, AspectRati 




- Graphics - 

torn = Modulator [L] [RandomBitSeq, Samplinglnterval T/10, NumberOf Curves -> i # 
ModulatingPulse -* FiltPulse] ; 

ListPlot [Im[ torn] , PlotJoined -> True] 




- Graphics - 
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ListPlot [Re [torn] , Plot Joined -> True] 




- Graphics - 

List Plot [Abs [torn] , Plot Joined -> True] 




- Graphics - 

xll = Table [TestPulse [6] [0] [t T] , {t, 0, 7 - 1/200, 1/200}]; 
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((Drop[xll, -1] -Drop[xll, 1]) 200 ) *2 / 200 // Apply [Plus, #]& 

0.92632 



TestPulse[6] [0] [t_] = LaurentC [4] [0] [t 5/7] 

3 5 t t t „, _ r , r. 3 5 t 



812500 



1. Sin[i/r[4, 

sinW 4, -LL]] 



5 t 



812500 



It := 5; 

PhaseFunction[t_] = 0L,t; 

phasepoints = Table [{t, PhaseFunction[t T] } // N, {t, 0, L , 1/40}]; 
ListPlot [phasepoints. Plot Joined -> True] 

NIntegrate: rncvb : NIntegrate failed to converge to prescribed accuracy after 7 



recursive bisections in tl near tl 



-1. 48053x10' 



Nlntegrate: :ncvb : NIntegrate failed to converge to prescribed accuracy after 7 
recursive bisections in tl near tl = - 1 . 12255 x 10~ 6 . 

NIntegrate: rncvb :- NIntegrate failed to converge to prescribed accuracy after 7 
recursive bisections in tl near tl = - 1 . 08808 x 10~ 6 . 

General: : stop : 

Further output of NIntegrate: rncvb will be suppressed during this calculation. 

NIntegrate: : slwcon : Numerical integration converging too slowly; suspect singularity, value 
of the integration is 0, oscillatory integrand, or insufficient WorkingPrecision. 
If your integrand is oscillatory try using the option Method- >Oscillatory in NIntegrate. 




TestPulse2 [6] [0] [t_] = LaurentC [4] [0][t 6/7] 

3 '6t^n r,r. 3 6 t 



Sin [*[4, 

torn = Modulator [L] [RandomBitSeq, Samplinglnterval T/40, NumberOf Curves -> 1, 
ModulatingPulse -► Filt Pulse] ; 

Save [ "ModulatorData.m" , {T, L, RandomBitSeq, FiltPulse, torn}] 



\ 
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ListPlot [Im[tom] , Plot Joined -> True] 




000 



- Graphics - 

ListPlot [Re [torn]. Plot Joined -> True] 



000 




- Graphics 
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BeginPackage [ "Iiaur en t Functions 11 " ] 
T::usage = "This is the syrcbol period" 
BT: : usage = "This is the usual product" 
h::usage = "This is the raw gaussian pulse" 

: usage = "This denotes modulation index Pi" 
^::usage ="^[I/,t} is Laurents function" 

hFiltered: :usage = "This is the filtered gaussian pulse" 

PhaseAngle: :usage = "This function takes some time to calculate. The following 

code will draw a graph of the f unctionPhaseFunction[t_] =N [ £ L t ] ; 
phasepoints =Table [ {t , PhaseFunction[t T] } //N, {t , 0, L ,1/40}]; 
ListPlot [phasepoints, Plot Joined -> True]" 

S::usage="S = Sin[£] " 

J: : usage = "or = © J * ■ 

C: :usaga = "C = Cos [35] " 

M: :usage = "M = 2 L " 1 M 

Kodulationlndex: :usage = "Modulationlndex = h" 

LaurentS: :usage = "Laurents [L] [n] [t] = Ein[ ^[L,t + n T] ] / S" 

LaurentC: : usage = "LaurentC [L] [K] [t] is Laurents C K , t " 

AIphaKI: :use.ge = "AlphaKI [LL] [K,i] is Laurents a K ,i " 

LaurentLK: .: ur age = "LaurentLK[L] [K] gives the support of C K ,t " 

The start of Modulator Definitions 

ANKInitialStateSetUp: tusage = 
"ANKInitialStateSetUp [L] [K] [XnitBitSeq, AccumulatedPhase] sets up the sequence 
of prior states of A K/N that the modulator went thgough to get to the 
constelation point specified by AccumulatedPhase which is really Ao.o" 

AKN: : usage = "AKN[L] [K] [ {State, AccumulatedPhase} ] defines A K , N in terms of A 0 . N " 

ModulatingPulse: :usage = 
"The Pulse is assumed to have the following structure Pulse [L] [KJ [t] " 

NumberOf Curves: : usage = "The number of pulses used by the modulator" 

Samplinglnterval: : usage = 
"Samplinglnterval is the interval between samples of the output of the modulator" 

InitialState: :usage = " Initial St ate is the set of bits that are assumed 
to be present before i.e {a_ 1# a_ 2 , . . . } " 

StartingQuadrant: rusage = "StartingQuadrant = A 0 ,_iand is a number? 

Modulator: :usage = "Modulator [L] [BitSeq, Opts] assumes the following 
default options StartingQuadrant -+ 0, InitialState -* Table[l, {i, 
1,20}] , Samplinglnterval -+ T/32 , NumberOf Curves 4 , ModulatingPulse-*- 
LaurentC . The Pulse is assumed to have the following structure Pulse [L] [K] [t] " 

Start of the Receiver Functions 
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FiltPulse::usago= "The default pulse and is called by FiltPulse[L] [K] [t] ■ 

Sync Sajr.pl e: : usage = "Given that the campling interval is T/32,then sync 

sample has rang -16 to 16 and this moves the point used to demodulate " 

Receiver: : usage = n Receiver [L] [ Input Seq, Opts] assumes the following- 
default options {StartingQuadrant->0 / InitialState->{l # l / l / l,l / l / l / l / 
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1} , Samplinglnterval -> T/32 , KodulatingPulse -> 
FiltPulse, SyncSample - > 0, NumberOf Curves -> 2 " 

RaceiverProper: :usrige = "ReceiverProper [L] [ StartingQuadrant, XnitialState, 
SamplingInterval,ModulatingPulse, SyncSample, NumberOf Curves, InputSeq] 
is called by Receiver after all the options have been resolved" 

Begin [ n * Private * M ] 

Vl*og[2] 

a : = ■ 

2 7TBT 



h[t_] 



r 2 7T o T 



Ideally we would and did define the effect of the convolution of h[t] by the formula below. However, this version of 
Mathematlca gives an error. 



r r fT h[t - x] 
hfii [t„] := Release [Module [{x}, J <51r JJ 



-T h[t- x] 

hFiltered[PulseWidth_] [t_] := Module [ {xl, x2, x3), xl = 

, *r rr fT h[tl-r] , r PulseWidth PulseWidth T 
(H[#l, 4C]fc> [Table[{tl, J ^ - dx}, {tl , , _}]]. 

x2 = Interpolation [xl] ; x2 [t] ] 

£ : = N [Modulationlndex tt] 

C : = Cos [$] ; 
S := Sin[$] ; 
J := ( 0 // Chop) ; 
M : = 2 L ~ 1 ; 

PhaseAngle [L_] [t_] /;ts0:=0 

PhaseAngle[L_] [t_] /; t i LT i= 3> 

PhaseAngle [L_] [t_] : = 
PhaseAngle[L] [t_J = Module [ {xl, x2, x3, x4, x5, x6}, xl = hFiltered [3 L T] f tl - -^-1; 

f t2 T 
x2 = Table[{t2, ff J Evaluate [xl] dtl}, {t2, 0, LT, ^^~}] ; interpolation [x2] [t] ] 

^[L„, t_] /; 0< t < LT := PhaseAngle [L] [t] 
^[L_# t_] /; 2LT > t ^ LT := 5 r PhaseAngle [L] [t - I* T] 
We need to put this to avoid the function being extrapolated 

^ [L_ / t_] /; ! (0<t<LT) && i (2 L T > t i L T) : = 0 
LaurentS[L_] [n__] [t_] :=Sin[^[L, t + nT]]/ S 



IxLurentFunctions.nb 



3 



AlphaKI [LLJ [K_, i_] /; (0 < i < LL) && (0 <= K < 2 LI - 1 ) : = 

„ . r , - f KNum - x2 \ 
Module [{xl, x2, x2 , KNum}, xl := |^x2 = Mod [KNum, 2]; KNum = ; x2j ; 

KNum = K; x3 = Table [xl, {ii, 0, LL - 1} ] ; x3 [[ij ] 

LaurentLK [L_] [K_] : = 
Modulo [{xl}, xl = Table [L (2 - AlphaKI [L] [K, ii] ) - ii, {ii, 1, L - 1} ] ; Min[xl]] 

LaurentC[L_] [K_] [t_J /; 0 <= K < 2 L : = 

LaurentS[L] [0] [t] ]~J LaurentS[L] [ii + I, AlphaKI [L] [K, ii] ] [t] 

ii = l 

The start of the modulator function 

ANKInitialStateSetUp[L_] [K_] [InitBit Seq_, AccumulatedPhase.] : = 
Module [{xl, x2, x3, x4, x5, acuphase, initbitseq}, 
initbitseq = InitBitSeq; 
acuphase = AccumulatedPha.se; 
UpdateSeg : = 

Module[{}, xl = acuphase - Sum [initbitseq [ [i ] ] AlphaKI [I#] [K, i] , {i, 1, L -1}]; 
acuphase = acuphase - First [initbitseq] ; initbitseq = Rest [ initbitseq]; xl] ; 
Table [ UpdateSeq, {i, 1, Laurent LK [L] [K] }] ] 

AKN[L_] [K_] [{State_, AccumulatedPhaseJ] : = 
AccumulatedPhase - Sum [State [ [i + 1] ] AlpbaKI [L] [K, i] , {i, 1, L -1}] 

Options [Modulator] := {StartingQuadrant -* 0, InitialState -* Table [1, {i, 1, 20}], 
S amp 1 i ngl nt erva 1 -+ T/32, H-^aiberOf Curves -* 4, ModulatingPulse-* LaurentC} 

Modulator [L_] [BitSeq„, Opts ] : = 

Module [ 

{xl, x2, x3, x4, x5, x6, state, AccumalatedPhase, " seq, AKKState, Curves, Pulse}, 
xl = Sajaplinglnterval /. {Opts} /. Options [Modulator] ; 
- tate = InitialState / . {Opts} /• Options [Modulator] ; 
x3 = StartingQuadrant/. {Opts} /. Options [Modulator] ; 
x4 = Saxrplinglnterval / . {Opts} / . Options [Modulator] ; 
Pulse = ModulatingPulse /. {Opts} /. Options [Modulator] ; 
Curves = (NumberOf Curves / . {Opts} /. Options [Modulator] ) -1; 
AccumulatedPhasa = x3 ; 
seq = BitSeq; 
Table [ 

AKNState[K] = ANKInitialStateSetUp [L] [K] [state, AccumulatedPhase] , {K, 0, Curves}]; 
x5 := Module [{}, state = Join [ {First [seq] } , Drop [state, -1]]; 
AccumulatedPhase = AccumulatedPhase + First [seq] ; 
seq = Rest [seq]; 

Table [AKNState[K] = Join[{AKN[L] [K] [{state, AccumulatedPhase}]}, 

Drop [ AKNS tat e[K], -1]], {K, 0, Curves}]; 
x6[r_] = Sum [ Sum [ ( j)*™t.t. w [ t i* xu pui s © [K ] [r + i T] , 

{i, 0, LaurentLK[L] [K] -1}], {K, 0, Curves}]; 
Table [x6 [r] , {r, 0, T - x4, x4}]]; 
Table [x5, {kk , 1, Length [BitSeq] } ] //Flatten] 

Options [Receiver] := {StartingQuadrant 0, 

InitialState -v {1, 1, 1, 1, 1, 1, 1, l, l, l, l, i, l, i, i, i, i, lf lf 1} , Samplinglnter- 
ModulatingPulse -> FiltPulse, SyncSample -> 0, NumberOf Curves -> 2}; 

Receiver [L_] [InputSeq_, Opts ] : = 

Module [{xl, x2, x3, x4, x5, x6}, 

xl = StartingQuadrant/. {Opts} /. Options [Receiver] ; 
x2 = InitialState/. {Opts} /. Options [Receiver] ; 
x3 a Samplinglnterval /. {Opts} /. Options [Receiver] ; 
x4 = ModulatingPulse/. {Opts} /. Options [Receiver] ; 
x5 = SyncSample/. {Opts} /. Options [Receiver] ; 
x6 = NumberOf Curves / . {Opts} / . Options [Receiver] ; 
ReceiverProper[L] [xl, x2, x3, x4, x5, x6, InputSeq] ] 
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Rec eiver Proper [L_] [ StartingQuadrant_, InitialStste_, Samplinglnterval_, 
MooulatingPulse_, Sync S amp le_, NumberOfCurves_, Input Seer ] : = 

Module [{xl, sgn, ReceivedSeg, ExpectedValue, seq, ReceiveNext, D, 5T}, 
xl = T / Sampling-Interval; 

ReceivedSoq = Partition [InputSeq, xl] // Transpose //# f [Sync Sample + 1]]&; 
ExpectedValue = 

ModulatingPulse[L] [0] [ (Laurent LK[L] [0] / 2) T + Sync Samp 1 e Samplinglnterval] 
& = StartingQuadrant; 
sgn =0; 
10 = {}; 

soq = ReceivedSeq; 
ReceiveNext := 
Module [{xl, x2}, 

xl = ( (-l) BSn First [seq]) // Im; 
If[Abs[xl - E>rpectedValue ] <= Abs [xl + ExpectedValue ] , 

D = Join[ 3D, {I}]; = & + 1, D = Join[ E>, {-!}]; = -1 ]; 

seq = Rest [seq] ; sgn = Mod[sgn + 1, 2]]; 
Table [ReceiveNext, {i, 1, Length [ReceivedSeq] }] ; 
©I . . 

End [ ] 

EndPackage[] 



